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On the work distribution for the adiabatic compression of a dilute classical gas 
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We consider the adiabatic and quasi-static compression of a dilute classical gas, confined in a piston 
and initially equilibrated with a heat bath. We find that the work performed during this process is 
described statistically by a gamma distribution. We use this result to show that the model satisfies 
the non-equilibrium work and fluctuation theorems, but not the flucutation-dissipation relation. We 
discuss the rare but dominant realizations that contribute most to the exponential average of the 
work, and relate our results to potentially universal work distributions. 
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When a system is driven away an initial state of ther- 
mal equilibrium by a mechanical perturbation, the sta- 
tistical distribution of work for that process exhibits uni- 
versal properties. In particular, the exponential average 
of the nonequilibrium work is related to an equilibrium 
free energy difference P, Q , 



o o 
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Furthermore, the work distribution for such a process, 
and the corresponding reversed process, are related by 
the following work fluctuation theorem 0, Q : 



Pk{-W) 
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Here, W is the work performed during a given realiza- 
tion of the process 5] ; (3 is the inverse temperature of a 
thermal environment with which the system is initially 
equilibrated; A_F is the free energy difference between 
two equilibrium states, both at temperature (3~^ , cor- 
responding to the initial and final values of the exter- 
nal work parameter; p is the work probability distribu- 
tion; and the subscripts 'i^' and 'i?' distinguish conjugate 
forward and reverseprocesses, where necessary. (See 
Refs. H S S S H 11 for more details and Ref. for 
an overview of related entropy fluctuation relations.) If 
we subject the system to a cyclic process then AF — 
and Eq. [II reduces to a result derived by Bochkov and 
KuzovlevllQ. 

The discovery of these relations makes it interesting 
to find model systems for which the work distributions 
can be computed analytically [IS llil ll^ 111 [ij, [iSl 111 
Ir^ nil Exact results have recently been derived 

for the example of a piston moving at arbitrary speed 
against an ideal gas 0, Q, . Here we consider the 
somewhat different case of the quasi-static compression 
or expansion of a dilute (but not ideal) classical gas. This 
model was suggested in email correspondence to one of us 
(C.J.) by Prof. Seth Putterman, and has also appeared 



FIG. 1: A gas confined to a cylinder with a controllable piston 



in this setting in a preprint by Prof. Jaeyoung Sung |20j . 
Using elementary statistical mechanics, we derive a non- 
trivial but tractable expression for the work distribution 
piW), Eq.^2 and use this to verify and illustrate Eqs.Q] 
and [3 

Let us define the model more precisely. Consider the 
system shown in figure^ a dilute classical gas confined in 
a piston. We assume that quantum effects are negligible, 
that the particles do not have any important internal 
structure, and that they rarely collide with one another. 
Specifically, the mean free path between particle-particle 
collisions is finite (unlike in Refs. 0,0,01), but much 
greater than the characteristic distance between nearby 
particles. Initially, the piston is held fixed and the gas is 
brought to thermal equilibrium with an external, infinite 
heat bath. The bath is then removed, preventing the 
further exchange of heat across the walls of the container. 
The piston is then very slowly forced inward, performing 
work as it compresses the gas to a new volume. In the 
corresponding reverse process we start with the gas at 
thermal equilibrium with the final volume of the forward 
process and then we adiabatically expand the gas back 
to the initial volume. 

It is useful to define a reference process, during which 
the gas remains in contact with the reservoir, and thus at 
constant temperature, as it is compressed reversibly; IS.F 
m Eqs. Ellis the free energy change during this reference 
process. By contrast, during the adiabatic compression 
described above, there is a steady rise in the kinetic tem- 
perature of the gas. Thus, although the gas continually 
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FIG. 2: The work probability density, p{W), given by Eg. 1111 The solid lines are the work of compression and the dashed 
lines are the negative work of expansion. Note that the magnitude of the work is greater for compression than expansion. 
Each distribution obeys the nonequilibrium work theorem and each compression-expansion pair are related by the work 
fluctuation theorem @. A direct consequence of the latter, illustrated in the figure, is that the corresponding forward and 
negative reverse work distributions cross at W = AF il'21 . 



self-equilibrates due to particle-particle collisions, it is 
driven away from the isothermal sequence of equilibrium 
states defined by the reference process. 

As a first pass at this model, let us use simple ar- 
guments to verify Eq. ^ In three spatial dimensions, 
the average equilibrium internal energy of a dilute gas 
of N identical particles is E = 3iV/2/3, and the en- 
tropy is given by the Sackur- Tetrode equation, S/N ~ 
\n{V/NA'^) + 5/2. Here V is the volume of the box, and 
A = a/ f3h'^/2TTm is the thermal de Broglie wavelength {h 
is Planck's constant and m is the particle mass). The 
free energy F ^ E — f3^^ S is then 
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which satisfies the scaling law (for any cr > 0) 



(3) 



(4) 



When such a gas is prepared in thermal equilibrium, as 
described above, its energy E can be viewed as a random 
variable sampled from the canonical distribution. 



PiE;p,V) 
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Z{(3,V) 



g{E;V)e-^'^, 
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where g{E) is the density of states and Z — e"^^ is 
the partition function. Since the pressure of a dilute 
gas is p = 2E/3V, and its energy during an adiabatic 
process changes by increments dE = —pdV, it follows 
that the product VE^^"^ is conserved as we slowly change 
the volume from Vq to Vi The final energy is thus 

El ~ (Vo/Vi)^/^£'o, and the work performed is 



W = Et -En::^ aEn 



Vn 



2/d 
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where d = 3. Note that a (and therefore W), is pos- 
itive for compression and negative for expansion, and 
that — 1 < a < cxD. For expansion the negative work 
is bounded by the initial kinetic energy. 
Defining q = a + l = (Vo/Vi)'^/'^, we get 
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ZjqP, Vq) 
Z{f3, Vb) 



dEog{Eo;Vo)e~<^P''° 
qf3Fiqf3,Vo) - pF{l3,Va) 



But qF{qi3,Vo) = F{l3,Vi) (Eq.^jl, hence the right side 
is simply /3AF, as predicted by Eq. ^ 

Although the above analysis is simple, it has its draw- 
backs. Eq. 13 is a large- A'^ approximation, not an exact 
result. Moreover, we have invoked macroscopic, thermo- 
dynamic arguments in deriving Eq. |H| Such arguments 
are valid when the aim is to describe the typical behavior 
of a system, but become suspect in the present context, 
since the average {exp{—(3W)) is often dominated by real- 
izations during which the system behaves very atypically. 
Finally, to verify Eq.|2we must solve for p{W), which re- 
quires obtaining the density of states, g{E). Therefore, 
we now proceed with a more careful analysis. For the 
sake of generality, we allow the dimensionality of space 
to be an arbitrary integer d > I, rather than assuming 
d = 3. 

The density of states is the derivative of the function 
<f>(i?), defined as the number of energy states with energy 
less than E. For a dilute gas in the classical limit. 
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where k = dN/2, and r(fc) is the gamma function. On 
the right side of Eq. Q the first factor accounts for the 
quantum graininess of phase space, the middle factor 
/N\ counts the number of arrangements of N identi- 
cal particles in a volume V , and the last factor is the vol- 
ume of a dA^-dimensional hypersphere of radius yhnE 
in momentum space |2]| . Hence 



1 {2T:mf 



dE /i2fe ^! r(fc) 



(8) 



The partition function Z = j dE ge and free en- 
ergy F now follow by direct integration: 



F{p,V) = -^lnZ{f],V)=-^ln 
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(We recover Eq. O with the approximation InA^! 
N\nN - N.) Eqs.inilHl and together give us 



(9) 
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We now solve for W by invoking the quasi-static in- 
variance of ^(E: V) [22.1231. As discussed in similar con- 
texts in Refs. [3, Il4[ l24l l25l|. since the gas continually 
self-equilibrates, the value of $ remains constant during 
the process. From Eq.[7|we get $ cx (VE'^/^)^ , therefore 
El = {Vo/Viy/'^Eo, which again leads to Eq. El only now 
for arbitrary d > 1. Thus, using Eq. we get 



p{W) = J dEoP{Eo)S{W-aEo) 
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e --'-eiaW), (11) 



where the unit step function 9 guarantees that W has 
the same sign as a. We see that p{W) for adiabatic com- 
pression (positive a), and p{—W) for adiabatic expansion 
(negative a) , are gamma distributions with shape param- 
eter k = dN/2 and scale s = These distributions 
are illustrated in figure HJ Note that the work distribu- 
tion depends on the ratio of the initial and final volumes 
and not on the absolute volume. 

Our result for p{W) allows us to verify the fluctuation 
theorem, Eq. ^ Let ap = {Vo/Vif/'^ - 1 and an ^ 
(Vi/Vb)^/'' — 1 denote the values of a for the forward 
(Vb Vi) and reverse (Vi Vq) processes. Note that 
— QfR/ap — (Vi/Vb)^^'^, and a^^+a^^ — —1. Combining 
these identities with Eq. we obtain 
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where 

A^^ = F(/3,T^i)-F(/3,yo) - ^In^ = ^^ln(l + a) , 

p Vi fj 2 

(12) 

by Eqs. H This confirms Eq. |21 

The validity of Eq. ^ follows immediately from 
Eq. El 0], though it can also be verified by the direct 
evaluation of / p{W)e~^^ dW. An alternative approach 
is to use a cumulant expansion 0: 



(13) 



where ojj is the j'th cumulant of p{W). Using standard 
properties of the gamma distribution , we get 



dN 
~2~ 



hence 



ln(e-/'^) 



dN ^ {-ay 



dN 



(14) 



ln(l + a), (15) 



again confirming Eq. for this model (see Eq. 112(1 . 
Truncating this expansion after two terms yields 



Ai^ w {W) 



2^^ 



1 dN , 
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(16) 



where (W) = loi and — lOi are the mean and variance 
of piyV). Eq. ^|is just the fluctuation-dissipation rela- 
tion of linear response theory. Naively, we might expect 
this to be an excellent approximation for our model, for 
either of two reasons. First, we have assumed a quasi- 
static process, apparently keeping the system in the near- 
equilibrium regime where linear response theory ought to 
apply. Second, for iV ^ 1, the central limit theorem sug- 
gests a Gaussian distribution of work values, and for a 
Gaussian only the first two cumulants are non-zero. How- 
ever, Eq. El reveals that the truncated expansion Eg. 1161 
is vahd only when |a| <C 1, that is, when Vx « Vq. Why 
does the system not respond linearly for larger |a|? 

First, Ea. ll6l is valid for small excursions away from the 
reversible, isothermal reference process described earlier. 
In our adiabatic process, however, the kinetic tempera- 
ture of the gas changes substantially, hence the system 
strays far from the reference path, unless Vi « Vq (From 
Eq. El and Ei = aEo it follows that the system ends 
with a canonical distribution of energies at temperature 
(Vo/Vl)^^'^/?"^)- Second, while p{W) is nearly Gaussian 
in the region around its mean, the average of exp(— 
is dominated by work values deep in the lower tail of the 
distribution, where the central limit theorem does not 
apply. Thus we cannot invoke the central limit theorem 
to throw out the higher (j > 2) cumulants; indeed, the 
relative sizes of the cumulants are independent of N. 
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Our results also illustrate the Clausius inequality, 

W = ^^«>^^Ml + a) = Af^, (17) 

since a > ln(l + a) for all real a. 

Assuming many particles, the initial energy of the 
gas is almost always very near to the equilibrium av- 
erage energy, E = dN /2I3, implying that the work per- 
formed during a typical realization is W^^'^ « aE. [As 
a consistency check, Eq. ^2 verifies that the peak of 
p{W) occurs at aE + 0{N^)] However, the average 
of exp(— /3T4^) is dominated not by these typical real- 
izations, but rather by those for which the work falls 
near the peak of p(W^) exp(-/?W^) From Eq.[ni 

we find that this peak occurs at « ctE/q, where 

q = {VolViY/"^. Using Eq.El we conclude that the domi- 
nant realizations are characterized by energies 
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E. 



(18) 



These are realizations during which the system begins 
with a very atypical energy (E/q), but ends in a mi- 
crostate that is characteristic of thermal equilibrium at 
temperature P~^- This is a generic feature associated 
with the convergence of (exp(— /3VF)) |2^ . 

As a practical matter, it is often desirable to fit exper- 
imental data to an appropriate probability distribution. 
Where the work is smooth and unimodal the gamma dis- 
tribution may be a reasonable parametric choice, since 
it explicitly arises in this physical example, it obeys the 
appropriate symmetries, and it can model the skew typ- 
ical exhibited by work densities. Indeed in a recent ex- 
periment, it was found that the work density could be 
adequately fit to a Pearson type III distribution, which 
is a gamma distribution generalized with a location pa- 
rameter ^2^. However, the gamma distribution is clearly 
inadequate for this task in general, since the distribution 
has a sharp lower bound. It has long been speculated that 
the generalized Gumbel (aka generalized Fisher-Tippett) 
distribution may provide a universal work distribution, in 
the sense that the work density for many different process 
may limit towards this general form. 
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_^ ^ab(x-s)+ae'''^-''> 



r{a)b 



(19) 



Here, a, b and s are real parameters. This distribution 
with_positive integer 'a' arises in extreme value statis- 
tics 29] and has been proposed as a general model of 
non-equilibrium fluctuations [s^,!^- We originally spec- 
ulated that this may be a plausible general distribution 
for work densities since this is the maximum entropy dis- 
tribution, from the location-scale family, given the mean 
work and the mean Boltzmann weighted work (Eq. [T]l. 
Moreover, it has recently been shown that the general- 
ized Gumbel provides an excellent fit to work densities for 



various alchemical transformations |32| ■ It is therefore in- 
teresting to note that for small b the generalized Gumbel 
approximates the generalized gamma distribution. Con- 
sequentially, the results of this paper are compatible with 
the universal distribution hypothesis. 
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